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ABSTRACT 

We assemble and observe a sample of poor galaxy systems that is suitable 
for testing N-body simulations of hierarchical clustering (Navarro, Frenk, & 
White 1997; NFW) and other dynamical halo models (e.g., Hernquist 1990). 
We (1) determine the parameters of the density profile p{r) and the velocity 
dispersion profile a p (R), (2) separate emission-line galaxies from absorption-line 
galaxies, examining the model parameters and as a function of spectroscopic 
type, and (3) for the best-behaved subsample, constrain the velocity anisotropy 
parameter, f3, which determines the shapes of the galaxy orbits. 

Our sample consists of 20 systems, 12 of which have extended x-ray emission 
in the ROSAT All-Sky Survey. We measure the 877 optical spectra of galaxies 
brighter than mj pa 15.4 within 1.5/i _1 Mpc of the system centers (we take 
Hq = 100/i km s~ x Mpc -1 ). Thus we sample the system membership to a 
radius typically three times larger than other recent optical group surveys. The 
average system population is 30 galaxies, and the average line-of-sight velocity 
dispersion is ~ 300 km s _1 . 

The NFW universal profile and the Hernquist (1990) model both provide 
good descriptions of the spatial data. In most cases an isothermal sphere is ruled 
out. Systems with declining a p (R) are well-matched by theoretical profiles in 
which the star-forming galaxies have predominantly radial orbits (f3 > 0); many 
of these galaxies are probably falling in for the first time. There is significant 
evidence for spatial segregation of the spectroscopic classes regardless of a p {R). 
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1. Introduction 

Groups of galaxies have attracted astronomers since the days of Shapley (1933) and 
Zwicky (with Humason, 1960). Although the division between clusters and groups must 
have seemed natural early on — comparing the "rich" Coma Cluster, for example, with the 
"poor" Stephan's Quintet — today the distinction is hardly obvious. A host of group catalogs 
(e.g., Albert, White, & Morgan 1977; Hickson 1982; Huchra & Geller 1982; Morgan & 
Hartwick 1988; Ramella, Geller, & Huchra 1989; Gourgoulhon, Chamaraux, & Fouque 1992; 
Ramella, Pisani, & Geller 1997; Trasarti-Battistoni 1998) make it clear that "groups" are 
often as well-populated as Abell's (1958) clusters (to the same apparent magnitude limit), 
and are sometimes parts of richer, perhaps bound systems (e.g., Ramella et al. 1994; Barton, 
de Carvalho, & Geller 1998). The potpourri of adjectives with overlapping meanings still 
used to describe groups and clusters — "poor," "rich," "compact," and "loose" — is a sign 
that the relationship among systems of galaxies across the full extent of their spectrum is 
not yet well understood. 

Much of the literature on "groups" deals with an important initial difficulty in studying 
poor systems of galaxies: optically selected catalogs are likely to contain large numbers 
of apparent systems which are actually chance superpositions even though the member 
galaxies have similar redshifts (e.g., Walke & Mamon 1989; Hernquist, Katz, & Weinberg 
1995). Detecting thermal bremsstrahlung x-ray emission from gas in the group potential 
increases the probability that galaxies in an overdense region form a bound system, (e.g., 
Schwartz, Schwarz, J., & Tucker, 1980; Ponman et al. 1996); however, if the large scale 
matter distribution is filamentary, the x-ray emitting gas as well as the galaxies may 
be projected onto a deceptively overdense region (Ostriker, Lubin, & Hernquist 1995). 
In conjunction with x-ray observations, obvious remedies against confusing the chance 
alignments with bound systems are to sample the overdense regions more deeply for 
additional associated galaxies (e.g., Zabludoff & Mulchaey 1998), and to identify additional 
group members within the neighboring few Mpc of the overdense region (e.g., Ramella et 
al. 1995). 

There is a rich corpus of theoretical work on the dynamics of low-velocity dispersion 
systems. N-body simulations of isolated groups initialized in virialized states (Bode, 
Cohn, & Lugger 1993; Bode et al. 1994; Athanassoula, Makino, & Bosma 1997) indicate 
that dynamical friction plays an important role in the evolution of groups, and that the 
formation of a central, massive galaxy which continually accretes smaller group members 
is rapid. In this theoretical scenario dynamical equilibrium is out of the question because 
of the large and roughly constant frequency of galaxy interactions. The cosmological 
simulations of Nolthenius, Klypin, & Primack (1997) similarly show that only halos with 
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mass > 1O 14 M are guaranteed to be virialized after a Hubble time; dynamical equilibrium 
is a poor description of less massive systems. However, Bahcall, Gramann, & Cen (1994) 
use their N-body simulations to claim that the velocity distribution of galaxies in groups 
is nearly Maxwellian, and Frederic's (1995b) N-body code provides evidence that certain 
robust virial mass estimates applied to groups are accurate and unbiased. Whatever the 
truth of the matter, equilibrium models often form the basis for conclusions about the mass 
distribution in groups (e.g., Persic & Salucci 1992; Dell'Antonio, Geller, & Fabricant 1995; 
Pedersen, Yoshi, & Sommer-Larsen 1997; Zabludoff & Mulchaey 1998). 

Literature which links real groups with theoretical work is less abundant. There has 
been some comparison of N-body simulations with observed group luminosity functions 
(LFs); Moore, Frenk, & White (1993), for example, find good agreement between group LFs 
and the results of their N-body simulations. Nolthenius, Klypin, & Primack (1994,1997) 
use other catalog statistics, such as the fraction of galaxies identified as group members, to 
constrain the cosmological parameters from N-body simulations; but Diaferio et al. (1999), 
whose simulations distinguish galaxies from dark matter halos, are able to reproduce the 
LFs of group catalogs largely independently of the cosmological parameters. 

For lack of adequate group data, however, direct comparisons of real systems with 
theoretical density and velocity dispersion profiles exist only for rich clusters (e.g. Carlberg 
et al. 1997a,b). The goal of this paper is to model the internal dynamics of low- mass 
systems of galaxies by providing a catalog suited to the task. 

We compare the optical properties of poor groups of galaxies with the Navarro, Frenk 
& White (1997) universal profile, which fits cosmological simulations of dark matter halos 
over many orders of magnitude in halo mass. We also examine the Hernquist (1990) 
profile, first derived to fit elliptical galaxies. Throughout the analysis we divide the system 
populations according to spectroscopic type by classifying the individual spectra as either 
absorption- or emission-line dominated; in this way we model the kinematics of galaxies 
with predominantly older stellar populations separately from those which have active star 
formation. 

To test the models, it is necessary (1) to have a large enough (« 30) galaxy membership 
per system and (2) to sample the membership well outside the core in order to cover the 
tails of the theoretical distributions sufficiently Previous catalogs of poor systems fulfill 
either the first (Pildis et al. 1995; Zabludoff & Mulchaey 1998) or the second (Ramella et 
al. 1995) criterion, but not both. In §0 we describe our assembly and observations of a 
catalog of groups, each with membership complete to roughly the same apparent magnitude 
(itir = 15.4) within 1.5k" 1 Mpc of the system center; in we discuss an objective method 
for spectral classification of the member galaxies. 
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We identify system members in §0, and in §|5| we determine the line-of-sight velocity 
dispersion profiles. We examine the behavior of the basic dynamical properties as a function 
of spectroscopic type in §|| and fit theoretical models in §[7|; in §[5] we summarize. 
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2. The Data 

We assemble and analyze a sample of galaxy systems, the Deep Optical Catalog or 
"DOC" hereafter. Two larger surveys are the sources of the DOC: (1) the Ramella et al. 
(1999) catalog of galaxy systems selected optically from the combined CfA-SSRS2 redshifts 
surveys (the "CSOC" hereafter; see below) and (2) the ROSAT All Sky Survey — Center for 
Astrophysics Loose Systems (the "RASSCALS" hereafter; Mahdavi et al. 1999), an x-ray 
emitting subset of the CSOC. 

Figure [I] shows the relationship among the DOC, the CSOC, and the RASSCALS. All 
the systems in the DOC and the RASSCALS also belong to the CSOC; however the DOC 
contains more galaxies per system, and has a high quality optical spectrum for each galaxy. 

2.1. Optical Source Catalog 

The CfA-SSRS2 Optical Catalog (Ramella et al. 1999) is a sample of galaxy systems 
drawn objectively from two complete redshift surveys. It includes systems in a wide variety 
of dynamical states, from groups with ~ 5 members to the Coma cluster. 

The Center For Astrophysics Redshift Survey (Geller & Huchra 1989; Huchra et al. 
1990; Huchra, Geller, & Corwin 1995; Falco et al. 1999 CfA) and the smaller Southern Sky 
Redshift Survey (Da Costa et al. 1994; Da Costa et al. 1998), both complete to m B = 15.5 
(rriR ~ 14.4), served as sources for the CSOC. Systems in the CSOC bear the label of the 
redshift survey in which they are located; the names of systems we describe here begin 
with either NRG (northern CfA) or SRG (southern CfA). The portion of the CSOC we use 
covers more than one sixth of the sky in two separate sections: (l)NRG, a 2 ooo = 8-5 hours 
to 17 hours and 5 2 ooo = 8. 5° to 43.5°, (2) SRG, 

a 2ooo = 21.5 hours to 3 hours and #2000 = 

-2° to 32°. 

Ramella et al. (1999) use a two-parameter method to construct the CSOC. Huchra & 
Geller (1982) first described this friends-of- friends algorithm (FOFA) for use with redshift 
surveys, and Ramella, Pisani, & Geller (1997; RPG) revised and applied it to the NRG 
data. The FOFA is a three-dimensional algorithm which identifies regions with a galaxy 
overdensity Sp/p greater than some specified threshold. A second fiducial parameter, Vo, 
rejects galaxies in the overdense region which are too far removed in velocity space from 
their nearest neighbor. The N-body simulations of Frederic (1995a) show that the Huchra 
& Geller (1982) detection method misses few real systems, at the cost of including some 
spurious ones. Ramella et al. (1999) apply the FOFA to the combined NRG and SRG 
redshift surveys with Sp/p = 80. 
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The CSOC contains 256 systems with 3000 km s _1 < (v) < 12000 km s _1 and an 
average of 9 members galaxies per system. The low velocity cutoff rejects systems which 
cover a large area on the sky, and which may be affected by motions in the Local Group. 
The median recession velocity for systems in the CSOC is 7000 km s _1 ; therefore the effects 
of cosmology and evolution are negligible throughout the sample. 

2.2. The X-Ray Data 

Mahdavi et al. (1997; MBGR) describe a technique to search the ROSAT All-Sky 
Survey (Voges 1993; RASS) for low surface brightness x-ray emission from systems of 
galaxies. Mahdavi et al. (1999) apply the method, which uses optical galaxy positions to 
limit the sky area searched for emission, to all the groups in the CSOC with 5 or more 
members. They use a newly processed version of the survey (RASS II; Boiler et al. 1998), 
which corrects effects leading to a low detection rate in the original reduction (RASS I). 

The MBGR approach sums x-ray photons only within the union of all projected 0Ah~ l 
Mpc regions around each system member identified in the CSOC. This "position template" 
method has an advantage over traditional detection algorithms: it operates on unbinned, 
unsmoothed, raw counts, and as a result can detect much lower surface-brightness emission. 
The drawback of the MBGR method is that it can measure only the total integrated 
emission within the position template; any analysis of points sources like active galactic 
nuclei (AGN) must take place post facto with traditional techniques. 

Using the position template method, Mahdavi et al. (1999) find that 101 of the 
CfA-SSRS2 systems have statistically significant x-ray emission in the RASS. In each case, 
they evaluate the probability P po int that the x-ray emission is due to a point source such as 
an active galactic nucleus (AGN). P po i n t results from a Kolmogorov-Smirnov Test (see §0) 
comparing the projected radial profile of the emission with that of the ROSAT PSPC point 
spread function. Systems with a high likelihood of having extended x-ray emission due to 
bremsstrahlung of hot gas in a potential (P po int ~ 0.03) make up the catalog of ROSAT All 
Sky Survey — Center for Astrophysics Loose Systems (RASSCALS). 

2.3. The Current Sample of Groups 

The Deep Optical Catalog (Table [1|) consists of 20 systems. Of these, 12 are selected 
randomly from the subset of CSOC groups which (1) have extended (P po int < 0.003) x-ray 
emission, (2) have velocity dispersions o v < 700 km s _1 , and (3) are observable from the 1.5 



- 8- 



meter Tillinghast Reflector on Mt. Hopkins, Arizona. The remaining systems were selected 
randomly from the subset of CSOC groups which fulfill (2) and (3) above, but either have 
pointlike or no x-ray emission. The systems with extended x-ray emission have a higher 
likelihood of being bound than the rest, but the presence of this emission is no guarantee 
that the groups are virialized. 

Because the CSOC systems are identified with objective overdensity and velocity 
difference criteria, the DOC is a representative sample of groups without implicit 
morphological selection biases. Our sample includes two Abell clusters, A779 and A1185. 
Figures Q-|6| show the RASSCALS x-ray data for this sample. 

MBGR analyzed the x-ray data for 10 of the DOC systems. However, the RASSCALS 
x-ray data differs somewhat from that of MBGR because (1) in MBGR the fields were not 
large enough to capture all the x-ray emission from some of the systems at low redshift, and 
therefore MBGR did not detect some of the RASSCALS, and (2) the RASSCALS x-ray 
data is from the more current reduction of the ROSAT All Sky Survey (RASS II). Mahdavi 
et al. (1999) discuss these issues in greater detail. 

We used the high-throughput FAST spectrograph on the Tillinghast (Fabricant et al. 
1998), which, with a 300 line mm -1 grating and a 3" slit, has a resolution of ~ 1.5 A, and a 
spectral coverage of 3940 A, centered at 5500 A. For z < 0.05, the redshift neighborhood 
of interest for this work, the FAST allows us to measure the equivalent widths of emission 
lines which are important markers of star formation, from [O II] at 3727 A to Ha at 6563 A. 

For each system we obtain digitized Palomar Observatory Sky Survey (POSS) positions 
for all galaxies with vtlr ^ 15.4, which lie within a projected distance of l.5h~ l Mpc of the 
system center as defined in the CSOC. We use the FOCAS (Valdes et al. 1995) program 
to separate stars from galaxies, but we also check each field by eye to make sure FOCAS 
has not neglected any bright elliptical or edge-on spiral galaxies. We calibrate the FOCAS 
magnitudes to the R magnitude system by (1) finding all CSOC galaxies with measured B 
band magnitudes in each field, (2) applying a standard B — R = 1.1 mag color conversion, 
the average for elliptical and spirals (Frei & Gunn 1994) to each CSOC galaxy, and (3) 
using a linear fit to the corresponding FOCAS magnitudes to obtain the transformation for 
all galaxies in each field. We estimate that the scatter around the tur = 15.4 completeness 
limit in each system is ~ 0.3 mag, and that from system to system it is ~ 0.5 mag. 

After obtaining the positions, we measure the spectra of all galaxies with unknown 
redshifts, and remeasure the spectra of galaxies with a known CSOC redshift in the interval 
\cz — (v)\ < 4000 km s" 1 , where (v) is the recession velocity of the system under study. 
Therefore the members of all the systems in the DOC, and most of the foreground and 



- 9- 



background galaxies, have high-quality FAST spectra complete to ph 15.4. 

There are a total of 877 new galaxy spectra in the DOC; 154 of these are remeasurements 
of CSOC galaxies. We analyze the data with the spectroscopy package of the Image 
Reduction and Analysis Facility (Valdes 1992; IRAF), and derive redshifts using the RVSAO 
program (Kurtz et al. 1992; Kurtz & Mink 1998), which applies cross-correlation techniques 
to log- wavelength binned spectra (Tonry & Davis 1979). We use cross correlation for both 
emission and absorption spectra. 

Table ^ lists the velocities. A sample page is included here; the entire table is available 
in electronic form. 
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3. Spectroscopic Classification 

We use the optical spectrum of each galaxy to decide whether the galaxy is "absorption- 
dominated," that is, consists mostly of an older stellar population (with metallic absorption 
lines that dominate the optical spectrum), or whether it is "emission-dominated" (with Ha, 
[O II] and other lines indicative of active star formation). 

To classify each of the new FAST spectra, we first de-redshift it to the heliocentric rest 
frame z — 0. We fit a cubic spline to determine the contiuum emission, which we subtract 
from the spectrum. Using only the portion between 3600-6900 A, binned into 2200 channels 
of width 1.5 A each, we fit the linear combination C e T e + C a T a to the spectrum, where C e 
and C a are the dimensionless emission and absorption coefficients determined from the fit, 
and T e and T a are the template emission and absorption spectra, respectively. 

These templates (Figure |7|) are the same ones we use to determine the redshift with 
RVSAO cross-correlation methods. Kurtz & Mink (1998) describe the template construction 
procedure, which for T e consists of averaging the profiles of the emission lines identified in a 
large number of emission-dominated galaxy spectra, and for T a involves averaging a large 
number of pure absorption spectra. The template spectra have no continuum emission. 
Before fitting the templates to the FAST data we extract the 3600-6900 A region. We bin 
the templates in the same way as we bin the data. 

We use the generalized least squares procedure (Press et al. 1992) to fit the 738 galaxy 
spectra with 500 km s -1 < cz < 15000 km s _1 . The distribution of x 2 / u f° r the fits has a 
mean value of 1.49. Figure |] shows shows the distribution of the parameters C e and C a . 
The 68.3% confidence errors are derived from the fit covariance matrix. 

In Figure [8], the absorption-dominated galaxies cluster at the lower right corner, 
where C e is zero. There is a smooth transition from this high density region to the 
emission-dominated galaxies, which rise towards the upper left corner of the figure. The 
emission-dominated population shows no clustering, because absorption lines from the older 
stellar population are detectable in almost all the galaxy spectra; thus most of the galaxies 
have C a ~ 1. 

Separating the two populations of galaxies is not straightforward, because their 
distribution does not show a clear break in the transition zone. Suppose, however, that 
there is a curve on Figure [8] which divides the populations in some maximum-likelihood 
sense. We would like to find this curve as objectively as possible, keeping in mind that 
some arbitrariness is inevitable. Luckily, a few physical considerations constrain the nature 
of the division: 
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1. The curve must be a function of the absorption coefficient with positive or zero first 
derivative. For suppose two galaxies have the same C e , but different C a ; and that the 
one with the smaller C a is classified as absorption-dominated. Then necessarily the 
other galaxy must be classified the same way. A similar argument holds for C e . 

2. The curve must not classify galaxies with C e < as emission-dominated. 

3. The curve must not classify galaxies with "large enough" C e as absorption-dominated. 

We separate the two galaxy populations on Figure |8| with a straight line C e = aC a + b 
for simplicity. Condition (1) above requires a > 0. Condition (2), given our data, requires, 
b > —0.6a. Selecting C e = ^c e + ^ a c £ — 0.3 as "large enough" (where ^c e is the mean 
and ac e is the standard deviation of the 738 C e 's ), condition (3) requires a < 1.5. The 
permitted region for a and b is now small enough to allow efficient numerical exploration. 

To find the optimum values of a and 6, we iterate over the permitted values, each 
time classifying galaxies with C e < aC a + b as absorption-dominated, and the others as 
emission-dominated. We exclude two galaxies with anomalous coefficients C e ~ C a ~ 
from this procedure. For each pair of populations we compute the two-dimensional KS 
statistic d 2 (see §y), which is a measure of the degree of difference between two populations; 
the larger o?2, the more likely it is that the two populations have different distributions. The 
maximum value of d>2 therefore gives the maximum-likelihood combination of a and 6. 

Figure |9] shows a contour plot of d,2 in the region of maximum likelihood. The optimum 
values are a = 0.02 and b = 0.038; the fit parameters are stable (to within the uncertainties) 
as the sample size increases. Figure § shows the line which separates the two populations. 
About 53% of the galaxies lie below the line; these are "absorption-dominated" galaxies; 
the rest are "emission-dominated." 
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4. Basic Dynamical Properties 

Pisani (1993) developed a non-parametric, scale-independent method which picks out 
structure in a one- dimensional data set. Without making use of a histogram, it estimates 
the probability distribution function (PDF) underlying the data as a kernel fka which is a 
sum of Gaussians with widths optimized to the local sampling rate: 

1 N 

fka{v) = —2_,K{cZ h S h v) (1) 



K(czi,s h v) = -== exp 



2ns 



(V - CZj) 2 



si 



(2) 



Here N is the number of galaxies in the field, each with a redshift Zj. The kernel is 
"adaptive" because s, changes for each data point. 

The adaptive kernel fka is a maximum-likelihood estimator of the system PDF. Pisani 
(1993) uses numerical simulations to show that fka-, a function with essentially N degrees of 
freedom, can accurately reconstruct even the most asymmetric and non-Gaussian velocity 
distributions. The kernel fka is particularly useful, because it provides an estimate of the 
probability that each galaxy in the field is a member of the system. We can therefore 
remove interlopers by discarding all galaxies with a large enough probability of belonging to 
the foreground or the background. Furthermore, whenever f^ a detects multiple peaks in the 
velocity distribution of a particular system, they are statistically significant substructures 
in velocity space. 

Figures pT0|— pT2 show histograms of the field of each system alongside fk a for that field. 
Our procedure for determining system membership is: 



1. Select all galaxies within ±4000 km s 1 of the system velocity as listed in the CSOC. 

2. Compute fka for these galaxies, as well as the probability that the ith galaxy is an 
isolated foreground or background object, P ie0 = K(cZi, Sj, czi)/ [fk a { cz i) N). Exclude 
all galaxies with P ie0 > 0.5 from further analysis. Note that a galaxy with P ie0 < 0.5 
is not necessarily a system member; it is just "unisolated." While the Pj £ o > 0.5 
cutoff is arbitrary, it corresponds well to traditional rejection algorithms, for example 
3-a clipping. 

3. Divide fka into v to t peaks by using its minima as delimiters. The galaxies are therefore 
also divided into v to t bins; the ith galaxy belongs to the z/th peak (j £ i/). 
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4. From ff. a compute the probability that the ith galaxy belongs to peak v. 



1 



K^cz^Sj.cZi) 



(3) 



Nfi 



5. Initially define the system membership as the peak v max with the most galaxies. 
However, if neighboring peaks have galaxies with a non-negligible probability of 
belonging to the main peak {P^ Vmax > 0.003), include these peaks in the system 
membership. 

Figures show sky plots of the members of each system, along with foreground and 
background galaxies and x-ray emission contours. 

After establishing the membership, we compute centers for each system. If there is a 
single peak in the x-ray emission and it is associated with member galaxies, we choose the 
center of the x-ray peak as the system center. If, however, the system has multiply-peaked, 
irregular, or no x-ray emission, we use the average right ascension and declination of the 
optically identified members. In a few cases this process moves the circle of radius 1.5/i -1 
Mpc away from the CSOC center, excluding some galaxies with measured spectra from the 
region of interest, and therefore from the membership. 

Finally, we compute the system recession velocity (v) and the line-of-sight velocity 
dispersion a p . In §[D] we test various robust methods of deriving these quantities; we find 



that the mean velocity v and the standard deviation y Var(u) are as effective as the robust 
estimators. 

There is a rich qualitative variation in the relationship between the system membership 
and nature of the x-ray emission: 

1. Of the 12 systems with significant extended x-ray emission, six (SRGbll9, NRGb032, 
NRGb045, NRGb244, NRGb251, and SRGb009) are characterized by a single x-ray 
emitting region which is relatively round and symmetric. Of these systems, five have 
either one or two bright elliptical galaxies at the peak of the x-ray emission: NRGb032 
(two), NRGb244 (one), NRGb251 (one), SRGb009 (two), and SRGbll9 (two). These 
systems increased their membership to three times the original CSOC count or better. 
The system without a bright elliptical, NRGb045, increased its membership by a 
factor of 1.6. Thus there is some indication that the presence of a bright elliptical 
galaxy at the x-ray core correlates with richness. The differences in the population 
increase are not due to uncertainty in the magnitude limit. 
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2. The remaining systems with extended x-ray emission, SRGb062, NRGb025, NRGsll7, 
NRGb247, NRGs385, and SRGb016, have irregular or multiple x-ray emission 
contours; all these systems more than tripled their original CSOC membership after 
our deeper sampling. 

3. Of the 8 systems without statistically significant extended x-ray emission, one 
(NRGb043) more than quintupled its original CSOC membership; this system is 
probably bound, and may emit below the RASS detection threshold. The rest 
(NRGb007, NRGb004, NRGb057, NRGsl27, NRGsl56, NRGbl81, and NRGs317) 
less than doubled the original count. Because these seven systems lack extended x-ray 
emission, and did not change their membership significantly after our survey, they 
may be artifacts of the FOFA algorithm. However, the deeper sampling does reveal 
a central condensation in both NRGb004 and NRGs317; it is therefore still not out 
of the question that these systems might be real groups which have x-ray emission 
undetectable in the RASS. 
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5. Velocity Dispersion Profiles 

5.1. On the Virial Radius 

In the literature describing N-body simulations of systems of galaxies (e.g., Navarro, 
Frenk, & White 1997), it is customary to show the dynamical properties as a function of 
r2oo, the radius which contains an overdensity 200p cr i t (z) (Navarro, Frenk, & White 1997; 
NFW), where p C rit(z) = 3Hq(1 + 2;) 3 /(87rG) is the critical density of an Einstein-de Sitter 
universe at a redshift z. 

Unfortunately, r2oo is not straightforward to derive observationally. One method 
(Carlberg et al. 1997a,b) is to assume that the virial theorem holds: 

M = aa 2 p r v G-\ (4) 

Here r v is the virial radius, G is Newton's gravitational constant, and a is a constant which 
depends on the orbit distribution of the system. Carlberg et al. (1997a,b) use a = 3, 
appropriate for an isotropic distribution. The next step is to assume that the mass inside 
a radius r is proportional to r: M{r) oc r 
redshift z implies 

|7rr| 00 x 200p crit (z) 
M v 

100r 2 2 00 ^(l + ^) 3 

r 200 



The (usually neglected) ~ 20% uncertainty in o v translates into the same error in r 2 oo- 
This loss of accuracy is significant considering that practically the only source of error in 
the projected radius R is the uncertainty in the determination of the center, when R is 
measured in Mpc. 

The following are drawbacks in deriving the observational r2oo from the virial theorem: 

1. If the orbits are anisotropic, a = 3 is incorrect by up to a factor of two. 

2. If M{r) is not proportional to r, the observational r2oo is inappropriate. NFWs 
"universal" density profile, for example, predicts M{r) oc r 2 for small radii and 



Then equation with a = 3 for a system at 



£ (6) 

T 



10 H 

aJl + z)- 3 / 2 , , 

—r-h- 1 Mpc. 

577 km s 



(7) 
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M(r) oc ln[r/(r c e)] for large radii. Carlberg et al. (1997a,b) show that this profile 
adequately describes rich clusters of galaxies; below (see §[7|) we find that it also fits 
our groups well, whereas an M(r) oc r profile is ruled out in most cases. 

3. Finally, the virial theorem may not apply to some of the systems in the sample. 
Even systems with extended x-ray emission are not guaranteed to be virialized. 
Substructure generally bloats, and infall generally compresses the line-of-sight 
velocity dispersion with respect to its equilibrium value. For example, the value of 
r 2oo — 1-2/zT 1 Mpc for NRGsll7 is probably too large, because this system has three 
subcondensations (Mahdavi et al. 1996). Similarly, r 2 oo — 0.12^i _1 Mpc for NRGb045 
is unlikely to be a true description of the overdensity radius, because o p = 70 km 
s _1 , from which it is derived, is smaller than the internal velocity dispersion of an L* 
elliptical galaxy! 



5.2. The Shapes of the Profiles 



Figures |10|-|T2] give a qualitative idea of the velocity dispersion as a function of projected 
radius. We divide the data for each group into three bins, containing the galaxies removed 
0-0.5 h~ l Mpc, 0.5-1.0 h^ 1 Mpc, and 1.0-1.5/t -1 Mpc from the center, respectively, and 
compute the velocity dispersion for each of these bins. 

Figure [13] shows the velocity dispersion profile cr p (R) for each system with N > 10 
members. We compute cr p (R) by sorting the members in R, and then computing the 
velocity dispersion of a moving group of 9 galaxies, from the center out to a distance of 
1.5/i -1 Mpc. Every ninth point in the Figure [13] is statistically independent. We compute 
the errors by performing the bootstrap analysis described in §[FJ on each bin. 

The behavior of u p {R) suggests that our sample consists of systems in a wide variety 
of dynamical states. The line-of-sight velocity dispersion profile, for example, rises in 
certain cases, falls in other cases, and sometimes just varies irregularly. The systems which 
have the most nearly constant u p {R) — NRGs317, NRGb057, and NRGb004 — are also the 
systems which have the fewest members; in these cases we may not be sampling the cores 
or outer regions sufficiently to detect a variation. Interestingly, these three systems also 
have x-ray emission which has a large probability of being due entirely to a point source 
(Ppoint — 0.055), and therefore have a smaller likelihood of being bound than the systems 
with clearly extended x-ray emission. 

To examine the variation among the a p (R) quantitatively, we perform the KS2D test 
(see §|Cp for the 105 unique pairs of systems with iV > 15 members. Table [3] shows the 
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test results separately for R in units of Mpc and r 2 oo- In both cases, for 41% of the unique 
system pairs in our sample, we can rule out the null hypothesis that <J p (R) is drawn from 
the same distribution at better than the 99.0% confidence level. 

The data show, therefore, that our galaxy systems are not, as a rule, in similar 
dynamical states, because their velocity dispersion profiles often differ significantly. If we 
had followed Zabludoff & Mulchaey (1998) or Carlberg et al. (1997), for example, in pooling 
all our data, we would have averaged away the distinctive rising, falling, and irregularly 
varying individual profiles. 

On the one hand, a few of our systems have profiles that qualitatively match recent 
N-body simulations of relaxed systems of galaxies as well as the analytical models we 
describe below (see §0). In SRGb009, SRGb016, NRGb032, SRGb062, and NRGb247, a p (R) 
resembles the velocity dispersion profiles which Crone, Evrard, Sz Richstone (1994) compute 
for simulated clusters in various cosmologies. The characteristic rise which all the Crone et 
al. (1994) models exhibit within the first w 0.2r 200 is present in SRGb009 and NRGb247, 
and cr p (R) declines for all five systems out to « 2r 2 oo, following the Crone et al. (1994) 
models. We will treat these systems, which have consistent distributions according to the 
KS2D test, as a special subsample ("Sample-D"). The x-ray emission in these systems is 
almost unquestionably extended (-P p0 mt < 0.001), another indication that they consitute a 
subsample which lends itself to clean dynamical modeling. 

On the other hand, across our sample, <J P (R) does not exhibit the regularity which it 
should if all our systems were in dynamical equilibrium, or if their members had similar 
families of orbits. Our results are consistent with recent observational studies of clusters 
of galaxies, (e.g. Girardi et al. 1996), which show that even a sample of systems with 
<jp ~ 1000 km s _1 exhibit both rising and falling velocity dispersion profiles within 1.5/i _1 
Mpc (« 2/i~ 1 r 2 oo)- Substructure, anisotropies in the true three-dimensional velocity 
distribution, and differences in the dark matter mass profiles, together or individually, may 
explain the variations among a p (R) for different systems of galaxies. 

We suggest that in the systems with declining a p (R) in the central region, the galaxies 
have mainly radial orbits, which bloat the radial velocity dispersion near the core, and 
shrink it at the outer edges. If this is the case, the galaxies could not have made many 
orbits; with a typical velocity of ~200-400 km s _1 , they travel only 2h~ 1 -4h~ 1 Mpc in a 
Hubble time, sufficient to make just one crossing of the system diameter. It is therefore 
possible that a significant fraction of the galaxies are falling towards the system center 
for the first time. In §[7| we will examine this scenario more rigorously, and consider the 
line-of-sight velocity dispersion profiles as a function of spectroscopic type. 
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6. Spectroscopic Segregation 

Elliptical and SO galaxies in clusters tend to be more centrally condensed than spiral 
galaxies, an effect related to the morphology-density relation in the nearby universe (e.g., 
Dressier 1980; Postman & Geller 1984; Lahav & Saslaw 1992; Whitmore, Gilmore, & Jones 
1993; Anderson 1996). Some authors (e.g., Carlberg et al. 1997) separate populations by 
color and refer to "red" and "blue" galaxies instead. Sometimes, elliptical or "red" galaxies 
also have a lower velocity dispersion than the remaining system members (e. g., Stein 1996; 
Colless & Dunn 1996; Mohr et al. 1996). Together, these two effects may imply that many 
systems perhaps contain virialized cores of galaxies with older stellar populations, with 
the star-forming "blue" spirals being more recent, infalling additions to the system. Our 
spectral classes correspond only roughly to morphological or color classifications; however, 
we seek to establish whether the stellar populations of galaxies vary with distance from the 
system center in a way similar to the morphology or the color. 

We compute (v), a p , and the average distance from the system center, R, as a function 
of spectroscopic type. Hereafter we use a superscript a with each of these quantities for 
absorption-dominated galaxies only, a superscript m for emission-dominated galaxies only, 
and no superscript for each system as a whole. We exclude NRGsl27, a system without 
x-ray emission, from our tests, because it contains only one emission-dominated galaxy. 

We compare the distributions of emission-dominated and absorption-dominated 
galaxies with the Student's t Test, which evaluates the probability that two populations 
have the same mean, and the F Test, which evaluates the probability that two populations 
have the same variance. When we apply the tests to individual systems, we use the 
measured velocities czi in km s _1 and projected distances Ri in Mpc from the center. 
Because scaling the data by a constant does not affect the outcome of either the Student's 
t or the F Test, it is only necessary to normalize the velocities and radii when we compare 
the absorption-emission populations across several systems in our sample. 

Therefore, when we pool together systems which have similar dynamical properties, we 
consider the normalized velocity Vi/a p in place of cz%, and measure R4 in units of r 2 oo as 
well as Mpc. We consider several pooled samples: "ALL" refers to all the galaxies which 
are system members; "Sample-A" refers to all the galaxies except the members of NRGsll7 
(A1185), a cluster with significant substructure (Mahdavi et al. 1996) which contains 14% of 
all the member galaxies, and therefore might significantly bias the test results; "Sample-D" 
refers to the systems with declining velocity dispersion (§ (5.2|) ; and "Sample-I" refers to 
all systems not belonging to Sample-D (and excluding NRGsll7). Table |] shows the test 



outcomes; Figure 14 shows the velocity and radius distribution of the galaxies in each 



system as a function of spectroscopic type. 
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The tests unambiguously show that our data recover a "spectroscopic segregation" 
which resembles the morphology-density relation if we identify absorption-dominated 
galaxies with the elliptical or "red" members. With better than 99.9% confidence, R a is 
different from R m when we pool all system members together and measure R in Mpc. On 
the average, the absorption-dominated galaxies are 0.14/i -1 Mpc closer to the system center 
than the emission-dominated galaxies. On an individual basis, nearly all the systems have 
R a - R m < 0; the two exceptions, NRGb007 and NRGb025, have R a = R m within the 
uncertainties. When we measure R in units of r2oo, the significance of the effect for the ALL 
sample drops to 99.0%. The physical units (Mpc) preserve the spectroscopic segregation 
relation more cleanly; the difference in physical scale among our systems is therefore smaller 
than error introduced into r 2 oo through the uncertainty in <r p ; see equation (|). 

The velocity dispersions of the absorption- and emission-dominated populations, 
however, show no sign of segregation. NRGsll7 is the only system for which P(F) (the 
probability that two samples are drawn from distributions with the same variance) is small. 
The results of the all the velocity tests for the combined subsamples are negative. 
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7. Comparison with Theory 



Here we compare the data with spherically symmetric density models from N-body 



subsample, Sample-D, contains groups with statistically consistent velocity dispersion 
profiles. These groups also have similar global velocity dispersions a p = 327-466 km s _1 , 
for which the variation in r^Q ial is comparable to its uncertainty. While in this section we 
focus most of our analysis on Sample-D, we also consider separately (1) the subsample of 
high velocity dispersion (a p > 350 km s _1 ) systems, and (2) the subsample of low velocity 
dispersion (cr p < 350 km s -1 ) systems. 

We consider spherically symmetric density profiles of the form 



where r c is the core radius, M c is the mass within the core radius, and S(a) ensures 
Jq c 4irr 2 p(r)dr = M c . The model with a = 3 is due to Hernquist (1990), who originally 
formulated it as a description of elliptical galaxies. The a = 2 model is the universal 
profile of Navarro, Frenk, & White (1996; NFW), who use it as a fitting function to 
their simulations of collapsed halos in various cosmologies. The a = 2 model is a good 
fit to the NFW simulations regardless of the cosmological parameters. The models have 
normalizations 5(3) = 3/8 and 5(2) = 31n(4/e)/2, respectively. We also consider the 
singular isothermal sphere (SIS), which has M(r) oc r and p(r) oc r~ 2 . 

Carlberg et al. (1997b) also fit these models to their sample of 16 high- luminosity x-ray 
clusters. Here we compare the models for the first time to a sample consisting primarily of 
groups of galaxies. Our approach differs somewhat from that of Carlberg et al. (1997b): 
below we compute the three-dimensional radial velocity dispersion profile analytically from 
the Jeans equation, rather than adopting a fitting function with four free parameters. 

For the sake of coherence among the model parameters, our notation for p(r) is 
somewhat different from that of Hernquist (1990), NFW, and Carlberg et al. (1997a,b). 
Because the NFW model has a logarithmically divergent cumulative mass M(r) as r — > oo, 
NFW define the total halo mass as M 2 oo = 200p cr i t (47r/3)r| 00 , where r 2 oo = cr c is the 
radius which contains 200 times the critical density of the universe, and c is the NFW halo 
"concentration." NFW find that the concentration varies inversely with the mass of the 
halo. We instead use the mass inside r c , which is more simply defined for both the NFW 
and the Hernquist (1990) models, for the profile normalization. 



simulations and the theoretical literature. We have shown (§ |5.2| ) that one particular 




(9) 
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7.1. Surface Number Density Profiles 

To fit the models to the optical data, we assume that galaxies trace mass. Then the 
three-dimensional galaxy number density is 

n(r) = ^p(r), (10) 

where N c is the number of galaxies within the sphere of radius r c . Our best fit iV c has a 
dependence on the rriR « 15.4 magnitude limit; the best fit r c should be independent of 
this limit if the mass-to-light ratio is constant throughout each system. The quantity N c is 
a superior fitting parameter to the core density n c = 3 N c / (47rr^) , because fitting for n c and 
r c , as Carlberg et al. (1997a,b) effectively do, results in an unacceptably large correlation 
between the parameters. In our experience, including the directly observable parameters 
as explicitly as possible in the models to be fit results in a smaller correlation among the 
parameters. 

The projection of equation ([K]) for a = 3 is the surface number density profile, 

2N r , r/_ 52 



and for a = 2 it is 



(r 2 -\ 



Z 3 (R) = x2 2 + R 2 ) X(R) - 3 , (11) 



^(R) = ; ; ° c 1 - X(R)] , (12) 



7rln(4/e)r2 (& - l) - 
where R = R/r c is the projected radius in units of the core radius, and 

X(R) = ^=1 (13) 
\JR 2 -l 

Note that X{R) is always real; for computation we use sech _1 i? = isec~ l R when R < 1. 
Also note that £ 3 (1) = 8N c /(15nr 2 ) and S 2 (l) = A r c /[31n(4/e)7rr^]. For the singular 
isothermal sphere, the surface density has the simple behavior Egig oc R~ x . 

To fit models, we use 20 bins, with an approximately equal number of galaxies ribin 
per bin per fit, so that the fractional Poisson error (y/riun h 1 10) _1 in each bin is roughly 
constant. However, if this fractional error is greater than 15%, we decrease the number of 
bins and increase n^n until the fractional error is smaller than 15%. We conduct two sets 
of fits: (a) measuring R in units of r^ 13,1 , given by equation (Q) (in which case r c = 1/c is 
dimensionless and equal to the inverse of the concentration); and (b) measuring R in units 
of Mpc. The x 2 function is well-behaved and has, in each case, a unique minimum; we 
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compute it over a fine grid around the minimum, and assume that the large-count Gaussian 
limit of the Poisson distribution describes our errors when we show the 68.3% and 95.4% 
confidence contours. Figures |H>HHj show the results of the fits for Sample-D, which we 
carry out in the log R — log S plane. We list the parameters for fits to other subsamples in 
Table |. 

The following are the outstanding properties of the fits: 

1. Both the Hernquist (1990) and NFW models provide good descriptions of the data. 
The isothermal sphere with M(r) oc r is usually the worst fit, and is often ruled out. 

2. The fits are of much higher quality when we measure R in units of Mpc than when we 
use r^Q iaJ ) the x 2 / u is always smaller. We conclude that the uncertainties discussed 
above (§ |5.1j) make a large contribution to the r^ 18,1 , and this quantity is therefore a 
poor estimator of the true r 2 oo for poor systems of galaxies. Classifying systems by 
the behavior of u p {R) (§ |5.2| ) is probably a more effective way to identify dynamically 
similar groups. 

3. The NFW concentration parameter c for Sample-D (a p = 327-466 km s -1 ) is 4.3-9.1 
(one-dimensional 95.4% confidence interval); the corresponding value for rich clusters, 
as derived by Carlberg et al. (1997a), is 2.3-7.7. In their simulations NFW find that 
low-mass halos should have larger concentrations that high-mass ones; our results for 
Sample-D, taken together with those of Carlberg et al. (1997a,b), are consistent with 
that prediction. 

4. The best-fit core radii r c for the high- and low-cr p systems would seem to contradict 
the NFW picture. If the low-cr p systems are truly less massive than their high-cr p 
counterparts, then NFW predict q > Ch as well as r 2 oo,z < r 2oo,c- Substituting the 
definition of the concentration, we have r c j < r c ^, which is not favored by the 
combined fit parameters listed in Table (^). We are nevertheless still faced with 
the fact that the NFW profile fits both subsamples well. Perhaps the internal 
inconsistency of the high- and low-cx p subsamples causes r c to behave this way. In 
§ [5.2| we showed that while Sample-D consists of systems with well-behaved velocity 
dispersion profiles, most other systems have profiles which are inconsistent. Therefore 
the conclusions drawn from Sample-D are more robust than those drawn from the 
high- and low-<r p subsamples. 

5. The fits support the results of our statistical tests in For all our subsamples, the 
best-fit r c is smaller for the absorption- than for the emission- dominated galaxies; the 
two-dimensional 68.4% confidence contours never overlap. 



-23- 



7.2. Velocity Dispersion Profiles 

The three-dimensional radial velocity dispersion v% for a spherical, nonrotating system 
is given by the Jeans equation (Binney & Tremaine 1987): 



IA. R) + = (14) 

par v ' r r z 



where r = r/r c and (3 = 1 — Vq/v% is the velocity anisotropy parameter. We solve the Jeans 
equation for the Hernquist (1990) profile, where 

4M c f 2 . . 

(r + l) z 

The Hernquist (1990) density profile yields a much simpler expression for than does 
NFW model; because both models fit our spatial data well, we only consider the former. In 
solving the equation we use constant (3, rather than assuming some functional form (3(r) 
from N-body simulations as Carlberg et al. (1997b) do. Because we are not dealing with 
rich clusters in this work, we would like freedom from the assumptions that the (3 profiles 
from simulations of clusters necessitate, e.g., that the velocity dispersion is isotropic near 
the core (/3 = 0) and that j3 achieves a maximum somewhere beyond the virial radius (Cole 
& Lacey 1996; Diaferio 1999). We therefore seek to establish the average value of (3 for our 
systems. 



For constant (3 and the Hernquist (1990) profile, equation (|14]) has the solution 

— AGM c (f + l) 3 f r 2/3_1 
f 2 ?- 1 J (f + 1 



V " = ~ 7. rW~l / 7^-^ ^ (16) 



While there is no analytic form for the above integral when (3 is a real number, it is possible 
to represent the solution as the following series (see §[E]), 

71 = ^(r~) 2 , (17) 
r c 

Af 00 

^ - (5 - 2 , )(r - +lf g^-+ir. (») 

where Co = 1. Note that we have set the integration constant equal to zero to supress 
divergent solutions as f — > oo. Equations pl7|-|T9| imply lim r ^o^r = 00 when (3 > 0.5, 



lim r ^ ^ 2 = GM c /r c when (3 = 0.5, and lim^o^r = when (3 < 0.5. The series always 
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converges for r > 0, because /3 < 1; it converges rapidly when either (3 < 0.5 or f > 1, 
and slowly when both /? > 0.5 and f < 1. In computing the series we use double precision 
numbers to minimize roundoff error, and add terms until we achieve a fractional accuracy 
of 10~ 6 . For the worst convergence case (r ~ 0, /3 — 1) we check the above series formula 
with the analytic solution 

r c 6r (r + 1) 

For this worst-case scenario, the series and the analytic solution produce identical results 
to 3 decimal places for f = 0.001 (which in practice never enters into consideration), and to 
5 decimal places when f > 0.01. 

The projected velocity dispersion profile is (Binney & Tremaine 1987) 



<^=mmL r^kCT (21) 



which we evaluate numerically by reformulating as 



9 ~ AGM C N C /■!/*/ ~ 9 9 n v(l/t) 2 t 2 dt . . 



Ttrpl{R)h y ( t + 1)3^1 _^2 t 

We integrate the inverse-square-root singularity by using the extended midpoint rule 
techniques described in Press et al. (1992; §4.4), which we modify to work with double- 
precision numbers. Note that once we specify r c from fits to the spatial data, the only free 
parameters in the theoretical u p {R) are (3 and M c . 

The combined projected velocity dispersion profiles for the low- and high-a p subsamples 
are too irregular to be good fits to simple density models. We fit only the combined profile 
for the well-behaved set of systems with declining a p (R) (Sample-D, described in §g). We 
fix r c at the value obtained from the spatial data, and fit (3 and GM C / (r c <Tp) as our two 
free parameters. The total mass of the Hernquist (1990) model results directly from our fit; 
M tot = 4M C . 



Figures 18-19 show the results of the fits, which have x 2 / u < 1 f° r both the emission- 
and absorption-dominated samples taken either individually or together. Therefore, a 
Hernquist (1990) profile, with a core radius taken from spatial data, fits the velocity data 
well, suggesting that Sample-D is well described by a spherically symmetric model with a 
simple velocity distribution function. 

It is perhaps more remarkable that the systems which have non-falling cr p (R) are still 
well fit by the simple Hernquist (1990) and NFW surface density distributions. Our data 
suggest that while a spherically symmetric mass distribution provides a good representation 
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of all our systems, the orbital families are irregularly populated, and (3(r,9,(f)) has a 
complicated behavior; it may not be represented as a spherically symmetric function even 
though the mass density is. This, in turn, is further evidence that our systems, with the 
exclusion of Sample-D, may not be in dynamical equilibrium. 

Because the best-fit (3 > for both absorption- and emission-dominated members 
in Sample-D, both sets of galaxies are probably dominated by radial orbits. The one- 
dimensional 95.4% confidence contour for the emission-dominated galaxies never extends 
below (3 = 0.2; we can therefore rule out orbital isotropy for them at that level. On the 
other hand, the orbits of the absorption-dominated galaxies are consistent with isotropy at 
the same level. 

Finally, we compute the masses traced by the galaxies. For Sample-D as a whole, the 
68.3% confidence interval on the total mass, not considering systematic effects, is 

M» = (1.2 ± 0.2) ft -xl0»M 9 ( M ^ :FT ) 2 . (23) 
The separate fits to the absorption- and emission-dominated galaxies yield 

M« t = (1.0 ± 0.2)*-' X 10"^( sj ^ Fr )*; (24) 

m S = (^•■r'xio'X^); (25) 

M m 

— ^ = 2.0 ±0.6. (26) 

lvl tot 

The absorption-dominated galaxies trace a mass that is entirely consistent with Sample-D as 
a whole. The emission-dominated members apparently trace a larger mass. The discrepancy 
in the mass estimates may be due to fitting a p (R) for a combined sample rather than on 
a case-by-case basis. It is also likely that equilibrium models do not provide an entirely 
appropriate description of the emission-dominated members. If the star-forming galaxies 
are on predominantly radial orbits, and are falling in for the first time as suggested by our 
timing arguments (§ |5.2|) , they probably have not reached dynamical equilibrium, and are 
not as well described by the Jeans equations as the absorption-dominated members. 
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8. Conclusions 



We analyze a sample of 20 galaxy systems drawn from a larger, objectively selected 
catalog of systems in the Center for Astrophysics redshift survey. Each system has a 
membership complete to tur ~ 15.4 out to a projected radius 1.5k' 1 Mpc from its center. 



1. Most of the systems have significantly more members after our observations than 
they had in the original source catalog (Ramella et al. 1997). Of the 12 systems with 
statistically significant, extended x-ray emission, 11 more than tripled their original 
membership. We conclude that our x-ray emitting systems have a high likelihood of 
being real, bound groups or clusters. Of the 8 objects without significant extended 
emission, one (NRGb043) quintupled its membership, and may have x-ray emission 
below the detection threshold of the ROSAT All Sky Survey; some of the other 
systems are unlikely to be bound configurations. 

2. We use fits to template optical spectra, along with the two-dimensional KS Test 
(§0), to classify the galaxies we observe into absorption- and emission-line dominated 
populations. We find evidence of spectroscopic segregation: the members with 
active star formation have consistently larger mean distances from the system center 
than the members with older stellar populations. Most systems have many more 
absorption- than emission-dominated galaxies in the projected regions encompassed 
by the x-ray emission contours. 

3. The line-of-sight velocity dispersion as a function of projected radius, cr p (R), is a 
useful tool for testing whether our systems are in similar dynamical states. A quarter 
of the sample exhibits declining velocity dispersion profiles similar to the those 
implied by N-body simulations (Crone et al. 1994; Navarro, Frenk, & White 1997). 
The remaining 75% of the systems exhibit flat, rising and irregularly varying profiles. 
The null hypothesis that a p {R) is drawn from the same distribution is rejected with 
better than 99% significance for 41% of all unique pairs of systems. 

4. Both the Hernquist model (1990) and the Navarro, Frenk, & White (1997; NFW) 
universal profile provide good fits to the spatial data, assuming that galaxies in our 
sample trace the system mass uniformly. The isothermal sphere is usually ruled out. 
Our best-behaved subsample (with intermediate velocity dispersions, a p = 327-466 
km s -1 ) yields a NFW concentration c = 4.3-9.1 (95.4% confidence interval). This is 
slightly larger than the range for rich clusters (Carlberg et al. 1997a,b), 2.3-7.7, and 
therefore consistent with NFWs prediction that the concentration should decrease as 
the mass increases. 
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5. We solve the Jeans equation for the Hernquist (1990) profile using a constant velocity 
anisotropy parameter /3; we project the solution numerically, and fit the theoretical 
cr p (R) to the best-behaved subsample. The theoretical profile provides a good fit, and 
the one-dimensional 95.4% confidence interval for (3 allows only predominantly radial 
orbits (f3 > 0) for galaxies with active star formation. Galaxies with older stellar 
populations have orbits which are consistent with isotropy to within the errors. The 
star-forming galaxies with predominantly radial orbits, especially those at a distance 
of ps 1.5/i _1 Mpc, probably have not made many crossings of the system diameter, 
since they move only 2h~ 1 -4h~ 1 Mpc in a Hubble time. Many of them may be falling 
in for the first time. 

6. Remarkably, the systems with irregular or rising velocity dispersion profiles have 
a spatial structure which is well-described by the Navarro, Frenk, & White (1997) 
and the Hernquist (1990) density profiles. Because these systems have surface 
number densities which are well-described by theoretical collapsed-halo profiles, but 
have velocity structure which no simple model can fit, they are probably bound 
configurations which are not in dynamical equilibrium. 

In future articles we plan a detailed spectroscopic analysis of individual system 
members to measure star formation rates, find galaxies which show evidence of truncated 
star formation, and model the relationship between the spectroscopic properties and the 
shape of the x-ray emitting region in our systems. 
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Appendices 
A. Definitions 

We use "system" to describe any set of galaxies which an objective structure- finding 
method places together as a unit. A system is "bound" if it is self-gravitating; it is 
an "artifact" if the objective method mistakenly placed physically unrelated galaxies 
together. The "membership" of a system is the set of all galaxies that have negative total 
gravitational energy in the center of mass frame; all other galaxies assigned to the system 
are "interlopers." 

A system with N members, each with a measured redshift z% and accompanying 
Gaussian uncertainty Azj, has a recession velocity (v), traditionally estimated by 
v = J2iC z i/N, where c is the speed of light. The ith member has a peculiar line-of-sight 
velocity V{ = [cz{ — (v)) / (1 + (v)/c). The global projected velocity dispersion of the system 
is a p , where traditionally is estimated by Var(t>) = J2i ( v i ~ ( v )) 2 /(N — 1). The standard 
deviation is defined as JVbx(v). In §|D| we consider robust estimators of (v) and <j v . 

We indicate three-dimensional distances from the system center with a lowercase 
r, while reserving the uppercase R for projected distance on the sky. In our article the 
"field" of a system is a circle on the sky with the same center as the system, and with 
radius 9 = 150 km s~ 1 /(w). In other words, the field is always large enough to include a 
R = 1.5/j. - 1 Mpc region around the system center. 



B. On Bootstraps With Measurement Errors 

When the probability distribution underlying the data, D(x), is a priori unknown — as 
is the case with the velocity distribution of system members — it is not possible to compute 
confidence intervals analytically for some parameter, a, derived from the data. The 
bootstrap is a numerical method for obtaining these confidence intervals by using the data 
set itself as an estimate of the underlying distribution (see, e.g., Lupton 1993). Classically, 
for a data set Xi of length N, 

1 N 

D(x) = -Y l 5(x-x i ) (Bl) 

i 

Here 5(x) is the Dirac delta function. To compute the confidence interval on a, one usually 
draws a large number of data sets of length N from D(x), and by recomputing a each time, 
derives its probability distribution. For the case that D(x) is a sum of delta functions, this 
procedure is tantamount to selecting points from the data set with replacement. 
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Unfortunately, by representing each data point as a delta function, D(x) neglects 
the uncertainty accompanying that measurement. The probability distribution of each 
individual measurement is not S(x — Xi), but some error distribution E(x) centered on the 
measurement. In the case of our redshift measurements, E(z) is a Gaussian with a standard 
deviation Azi « 10~ 4 . Therefore the correct bootstrap data probability distribution is 



Drawing deviates distributed as D(z) is a bit more tricky than selecting with replacement, 
as is possible with the delta function distribution. We use a rejection method, as outlined 
for example in Press et al. (1992, §7.3), with a constant comparison function over the range 
of the data. All the bootstrap confidence intervals for the parameters we derive therefore 
take the measurement uncertainties into account properly. 

Note that D(x), D(z) and the adaptive kernel estimator f ka , defined in equation (JJ), 
are only superficially related. The first two are minimal representations of the distribution 
of the data, constructed for the purpose of deriving alternate data sets drawn from a similar 
parent population. The adaptive kernel, on the other hand, attempts to overcome the 
effects of sampling as it builds its guess at the true shape of the parent distribution. 



The Kolmogorov-Smirnov Test (e.g., Press et al. 1992, §14.3) evaluates the null 
hypothesis that a pair of data sets are drawn from the same distribution. In the simple case 
that each data set is one dimensional, the correponding KS Test (KS1D) computes d\, the 
maximum deviation of the two data cumulative distributions. The a priori distribution of 
d\ itself is known, and therefore so is the probability of the null hypothesis. The KS1D test 
is particularly useful, because its estimate of the probability is independent of the original 
distributions of the data, which are often not known. 

We use a version of the Kolmogorov-Smirnov test (KS2D; e.g. Press et al. 1992, §14.7) 
which is applicable to two-dimensional data sets. Because the notion of a cumulative 
distribution is poorly defined for a two-dimensional data set, the corresponding deviation 
statistic, di-, is calculated somewhat differently. Suppose we are comparing two data sets, 
(x^ yi), i — 1 . . . N\ and (Xj, Yj), j = 1 . . . N 2 . For each (xj, yj) compute 





C. 



On the Two-Dimensional KS Test 



T>i = max 
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E 



x<Xi, y>y, 



Ni 



E 

X<Xi, Y>y, 



1 
N~ 2 



E 



x>x it y>y, 



Ni 



E 

X>x u Y>y t 



1 



(CI) 



Define T>~ similarly, replacing i with j, with X,-, and and K in the above equation. 



Then 



do = 



( N±N 2 
\Ni + N 2 



x 



max(T'.j) + max(Vj] 



(C2) 



To compute <i 2 we use the computer code which Press et al. (1992) implement on 
the basis of theoretical work by Fasano & Franceschini (1987). However, we do not use 
the Press et al. (1992) prescription to compute the probability of the null hypothesis, 
because it is only accurate for N ^ 20 data points. Rather we use the results of Fasano & 
Franceschini (1987) directly; they fit a third-order polynomial in three variables (sample 
size, sample correlation, and desired significance level) to their Monte Carlo simulations of 
the probability distribution of d 2 . They find that the fractional uncertainty of the resultant 
probabilities for the null hypothesis is ~ 5%. 

Note that the output dp of the Press et al. (1992) computer program is normalized 
differently from our notation for d 2 , which agrees with that of Fasano & Franceschini (1987): 



do = di 



f NxN 2 
\N 1 +N 2 



(C3) 



D. Comparison of the Robust and Classical Estimators of the Velocity 

Moments 

The classical estimators of the system recession velocity (v) and the line-of-sight 
velocity dispersion a p are the mean velocity v and the unbiased standard deviation estimator 



Var(t> ), respectively (see §0). However, Beers, Flynn, & Gebhardt (BFG; 1990) argue that 
the sample mean and the sample standard deviation are not, in general, efficient estimators 
of the properties of a parent distribution p(x). For example, if p(x) is not a Gaussian, x is 
in general a poorer estimator of (x), and Var(x) is a poorer estimator ((x — (x)) 2 ), than 
certain robust estimators which we compute alongside the classical estimators. 



In Figure EUI we show v alongside the sample median M, and the biweight location 
estimator, 

BiMean = M + (<g| - M) " ~ »' f ; (Dl) 
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^ c x MAD ' 1 j 

MAD = median (\ Xi - M\) . (D3) 

Here c is an arbitrary tuning parameter. Similarly, we list the standard deviation alongside 
the "Gapper," which uses the weighted gaps among the data points to compute the scale 
of the velocity distribution (see BFG for its definition); we also show the biweight scale 
estimator, 

rvq- 1/2 [E N<1 (^-M)(1-^) 4 ] 1/2 

BiSigma = n ' D4 

E| W |<1 (1 - M?) (1 - 5/xf 

Interestingly, these various estimators differ very little (Figure |20p. For all systems, the 
robust biweight estimates of the recession velocity and velocity dispersion lie well within the 
68.3% confidence interval of the classical estimators, the mean and the standard deviation. 
The median, which is not efficient in the presence of substructure (BFG), lies outside the 
confidence interval of the mean only twice. The even closer agreement of the estimators of 
cr p , a quantity somewhat more important for kinematic analysis than (v), is also noteworthy. 

In considering which of the two kinds of estimators, classical or robust, we should 
use, we note that the latter make use of a dimensionless tuning constant c, which BFG 
set to 6.0 for the biweight location estimator, and to 9.0 for the scale estimator. The 
biweights also require the subjective choice of an auxiliary estimator of scale; BFG use the 
mean absolute deviation ("MAD"). On the other hand the classical estimators contain no 
arbitrary parameters. 



We therefore adopt (v) = v and a p = \JYbi{v). All the quoted errors on the velocity 
moments are 68.3% bias-corrected bootstrap confidence intervals as described in §[B[ 



E. Integration of Constant-/? Jeans Equation 



Equation (jig ) is integrable through the identity, 



/ x m {a + bx) n dx 



x 



m+l , 



a + bx) r 



tin 



m + n + 1 



m + n + 1 



x m (a + bx) n - 1 dx. 



(El) 



Substituting a = b = 1, x 



2/5 — 1, and n = —5, and applying the above 



expression to the integral in equation [LJ| we obtain 

I { 5 



^2/3-1 
if + l) 5 



dr 



2(3-5 1 (f + l) 5 2/5-6 



6 



[f+lf 2/5-7 \(r + l) 7 



GE2) 
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(f+l) 5 5-2/?\ 6-2/3 



1 6/1 

(r + 1) + 7-2/3 \ (f+iy 



From the above relation, equations 17-19 follow easily. 
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Fig. 1. — Venn diagrams showing the relationship among the Deep Optical Catalog (This 
paper), the RASSCATS (Mahdavi et al. 1999), and the CTA-SSRS2 Optical Catalog (Ramella 
et al. 1999). In (a), the circles represent the set of all galaxy systems in each catalog. The 
area of each circle is proportional to the number of systems it contains; the areas of the 
intersections accurately represent the degree of overlap among the catalogs. In (b), the 
circles represent the set of all galaxies classified as system members in each catalog. The 
area of each circle is proportional to the total number of galaxies classified as members, and 
the areas of the intersections accurately reflect the number of galaxies the catalogs share. 
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Fig. 2. — Galaxy positions and x-ray contours. The x-ray emission contours begin at 1.33 
standard deviations above the level of the background, and increase by a factor of 1.5; dotted 
contours lie less than 3.0 standard deviations above the background. Galaxies with cz > 
15000 km s _1 are simply shown as crosses; they are never system members. Galaxies with 
500 km s _1 < cz < 15000 km s _1 are marked by circles if their spectra are absorption 
dominated, or by triangles if the spectra are emission-dominated. Each triangle or circle is 
filled-in if the galaxy is a system member, empty if it is not. 
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g. 3. — See caption for Figure [2|. 
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4. — See caption for Figure 0. 
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\. 5. — See caption for Figure ||. 
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Fig. 7. — The absorption and emission spectral templates, with a list of the lines present. 
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Fig. 8. — C a vs. C e for the 738 galaxies with 500 km s^ 1 < cz < 15000 km s -1 . Triangles and 
circles represent emission-dominated and absorption-dominated galaxies, respectively. The 
errorbars show the typical uncertainty in the classification; the straight line is the maximum- 
likelihood dividing line described in the text. The top right box is a magnification of the 
absorption-dominated region. 
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Fig. 9. — Maximization of the KS2D statistic (§0) for spectroscopic classification. Shown 
are contours of uniform d.2 = 11,12,12.6,13, and 13.4 in the a-b (slope-intercept) plane. 
Tickmarks indicate the downhill direction, and the solid thick line indicates the b > —0.6a 
boundary. The maximized value of d,2 is 13.44; a "+"marks its position. 
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Fig. 10. — The distribution of vt as defined in section §0. Histograms are coded (1) by 
line type: double-lined histograms show system members, heavy-lined histograms show 
foreground and background galaxies inside 1.5 h^ 1 Mpc of the system center, and dotted 
histograms show galaxies that are no longer inside this region after we recenter the system; (2) 
by shading: completely filled histograms show all galaxies within 0.5/iT 1 Mpc of the system 
center, cross-hatched histograms those within 0.5/i _1 -1.0/i _1 Mpc, and empty histograms 
those within 1.0/i -1 — 1.5/i -1 Mpc. On the right is a plot of fk a a s defined in equation flU). 



-47- 




Fig. 11. — See caption for Figure 
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Fig. 12. — See caption for Figure 
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Fig. 13. — Velocity dispersion profiles. The error bars show the 68.3% bias-corrected 
bootstrap confidence intervals. The number in parenthesis indicates the global line-of-sight 
velocity dispersion a p in km s _1 . The dashed vertical lines show the 68.3% confidence interval 
for r 2 oo as defined in equation (l|). 
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Fig. 14. — Velocity-distance diagrams. Dotted lines represent the mean projected distance 
and mean velocity of the emission-dominated galaxies (triangles), with accompanying errors; 
solid lines represent these quantities for absorption-dominated galaxies (circles). The number 
next to each system name is the probability from the Student's t test, that the emission- 
dominated and absorption-dominated galaxies have the same mean projected distance from 
the system center. 
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0.1 1 0.1 1 0.1 1 

R (tr 1 Mpc) 

Fig. 15. — Surface number density profiles and best fit Hernquist (1990) NFW, and 
isothermal (SIS) profiles for Sample-D. The top panels show fits with R measured in units 
of r 2 oo; the bottom panels show fits with R in Mpc. From left to right we have (a) the fit to 
the combined emission- and absorption-dominated populations; (b) the fit to the absorption- 
dominated galaxies only; and (c) the fit to the emission-dominated galaxies only. 
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Fig. 16. — Contours of the 68.3% and 95.4% confidence regions for the fits in the top panel 
of figure [15. 
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Fig. 17. — Contours of the 68.3% and 95.4% confidence regions for the fits in the bottom 



panel of figure [15 
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Fig. 18. — Fits of the line-of-sight velocity dispersion profile <j p (R) to the Hernquist (1990) 
model with constant (3. Shown are (a) the fit to the combined emission- and absorption- 
dominated populations; (b) the two populations separately, where solid error bars represent 
the absorption- and dotted error bars represent the emission-dominated members; (c) the 
fit to the absorption-dominated galaxies only; (d) the fit to the emission-dominated galaxies 
only. 
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Fig. 19. — Confidence contours corresponding to the fits in figure [18. 
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Fig. 20. — On the left for each system in our sample, we plot a stack of the three estimators of 
location, along with their 68.3% confidence intervals; shown are the mean (top), the median 
(middle), and the biweight (bottom). We shift the values for each system so that the mean 
is zero. Similarly, on the right we plot for each system a stack consisting of the standard 
deviation (top), the Gapper (middle), and the biweight (bottom) estimators of scale, shifted 
so that the standard deviation is zero. 
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Table 1. The Deep Optical Catalog: Basic Data 



Group 


a 


8 


N 


N 


N 


(v) 




(Tp 






p 


ID 


J2000 


J2000 


csoc 


Total 


DOC 


km s 


i 


km s 


— l 


Comments 


point 


SRGb062 


00:18:22.5 


30:04:00 


13 


84 


45 


6819 + 


57 
59 


379+ 


40 
42 


X,XC 


< 0.001 


SRGbll9 a 


01:56:13.8 


05:35:12 


8 


54 


28 


5466+ 


71 
67 


3671 


43 
43 


X,XC 


< 0.001 


NRGb004 


08:38:07.3 


24:58:02 


9 


36 


19 


8559+. 


61 
61 


2741 


50 
54 


X,XC 


0.587 


NRGb007 


08:50:29.9 


36:29:13 


6 


39 


6 


7541± 


1!) 
20 


27± 


15 
12 


OC 




NRGb025 b - c 


09:13:37.3 


29:59:58 


5 


35 


31 


6735+. 


84 
79 
6!) 
69 


4581 


70 
61 


X,XC 


< 0.001 


NRGb032 b ' d 


09:19:46.9 


33:45:00 


5 


67 


47 


67731 


4661 


67 
66 


X,XC 


< 0.001 


NRGb043 b 


09:28:16.2 


29:58:08 


5 


30 


26 


7891+. 


48 
57 


2551 


33 
33 


OC 




NRGb045 b 


09:33:25.6 


34:02:52 


5 


26 


8 


8181+. 


27 
31 


70l 


19 
18 


x,xc 


< 0.001 


NRGb057 b 


09:42:23.2 


36:06:37 


5 


24 


13 


6766^ 


52 
46 


1681 


24 
30 


x,xc 


0.684 


NRGsll7 b '° 


11:10:42.9 


28:41:38 


14 


89 


84 


97991 


71 

76 


6631 


42 
45 


x,xc 


< 0.001 


NRGsl27 b 


11:21:34.2 


34:15:31 


8 


12 


10 


104851 


59 
71 


2061 


36 
38 


OC 




NRGsl56 b 


11:44:54.5 


33:15:17 


7 


35 


28 


96401 


65 
62 


32ll 


35 
42 


x,oc 


0.111/0.210 * 


NRGbl81 b 


12:07:35.5 


31:26:32 


6 


18 


9 


67971 


57 
62 


1761 


35 
39 


OC 




NRGb244 a 


13:23:57.9 


14:02:37 


7 


32 


19 


69481 


61 

63 


2691 


34 
36 


x,xc 


0.003 


NRGb247 


13:29:25.7 


11:45:21 


12 


66 


39 


685ll 


66 
60 


3951 


40 
41 


x,xc 


< 0.001 


NRGb251 


13:34:25.3 


34:41:25 


5 


48 


26 


73461 


49 
56 


2651 


31 

35 


x,xc 


< 0.001 


NRGs317 


14:47:05.3 


13:39:46 


11 


46 


18 


88731 


78 
78 


33 ll 


36 
41 


x,oc 


0.055 


NRGs385 


16:17:43.9 


34:58:00 


8 


61 


53 


93521 


75 
80 


5631 


37 

40 


x,xc 


< 0.001 


SRGb009 


22:14:46.0 


13:50:30 


8 


46 


36 


7744+ 


66 
64 


39ll 


55 
43 


x,xc 


< 0.001 


SRGb016 


22:58:45.9 


26:00:05 


6 


58 


43 


73951 


52 
54 


3271 


33 
34 


x,oc 


< 0.001 



Note. — Nqsoc is the number of galaxies classified as members in the CSOC; iVrotai is the total number 
of galaxies to to# w 15.4 in the field of the system; Adoc is the total number of galaxies classified as 
members in this work. P po int is the probability that the x-ray emission in the field is due to a single point 
source. Letters in the "comments" column denote the following. X: has x-ray emission. XC: Position quoted 
is the x-ray centroid. OC: Position quoted is the mean optical position. 

a Also studied by Zabludoff & Mulchaey (1998). 

b Also studied by Ramella et al. (1995). 

c Contains Hickson Compact Group 37 (Hickson 1982). 

d Abell 779 (Abcll 1958). 

c Abell 1185 (Abell 1958). 

f Probabilities for the south and the north x-ray emission peaks coincident with member galaxies, 
respectively. 



-58- 



Table 2. Galaxy Positions and Measured Velocities 



Galaxy ID 


«2000 


^2000 


cz l 


SRGb062.01 


00:11:57.0 


+29:29:08 


28038 


SRGb062.02 


00:12:11.8 


+29:19:09 


7728 


SRGb062.03 


00:12:17.5 


+29:52:17 


6921 


SRGb062.04 


00:12:28.4 


+29:32:38 


6862 


SRGb062.05 


00:12:38.3 


+30:06:08 


6791 


SRGb062.06 


00:12:45.0 


+29:22:15 


10419 


SRGb062.07 


00:13:12.7 


+31:08:43 


14465 


SRGb062.08 


00:13:45.1 


+30:11:40 


7109 


SRGb062.09 


00:13:55.9 


+28:45:47 


6908 


SRGb062.10 


00:13:57.2 


+30:52:55 


4782 


SRGb062.11 


00:14:01.9 


+29:25:57 


7055 


SRGb062.12 


00:14:13.7 


+28:52:36 


7333 


SRGb062.13 


00:14:55.2 


+31:06:12 


24282 


SRGb062.14 


00:15:07.6 


+28:52:29 


38145 


SRGb062.15 


00:15:15.6 


+29:21:44 


7046 


SRGb062.16 


00:15:22.2 


+29:39:42 


6790 


SRGb062.17 


00:15:23.1 


+30:43:14 


14292 


SRGb062.18 


00:15:24.7 


+29:38:48 


21460 


SRGb062.19 


00:15:28.8 


+30:43:23 


13890 


SRGb062.20 


00:15:43.9 


+29:39:58 


6642 


Note. — This is a 


, sample listing; the full table 


is available electronically. All velocities 


are in km s _1 ; the 



typical uncertainty is 40 km s 
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Table 3. Kolmogorov-Smirnov Tests for a(R) 



(01) (02) (03) (04) (05) (06) (07) (08) (09) (10) (11) (12) (13) (14) (15) 



Note. — The results for R measured in units of Mpc are above and to the right of the dividing diagonal, 
while those for R in units of r2oo are below and to the left. For each unique pair of systems, the null 
hypothesis that <j{R) is drawn from the same distribution is rejected (1) at less than the 95.4% confidence 
level if there is only one star; (2) at better than the 95.4% confidence level if there are two stars; (3) at 
better than the 99.0% confidence level if there are three stars; and (4) at better than the 99.9% confidence 
level if there are four stars. 



Table 4. Comparison of Absorption- and Emission-Dominated Populations 



Velocity Projected Distance 



Group /abs P(F) allal P(t) M!_M1 P(F) P{t) Ra _ R e 

h^ 1 Mpc 



SRGb062 D 





.689 


SRGbll9 





.571 


NRGb004 





.632 


NRGb007 





.333 


NRGM25 





.613 


NRGb032 D 





.766 


NRGb043 





.538 


NRGM45 





,625 


NRGb057 





.385 


NRGsll7 





.560 


NRGsl56 





.571 


NRGbl81 





.444 


NRGb244 





.158 


NRGb247 £> 





.538 


NRGb251 





.654 


NRGs317 





.667 


NRGs385 





.717 


SRGb009 D 





.667 


SRGb016 D 





.535 


ALL 





.602 


Sample- A 





.609 


Sample-D 





.646 


Sample-I 





.585 



0.129 1.49 ±0.35 

0.805 0.94 ±0.22 

0.914 0.99±0.55 

0.727 0.52 ±0.87 

0.264 1.39 ±0.46 

0.804 0.96 ±0.26 

0.775 1.09 ±0.32 

0.110 0.39±0.26 

0.636 1.20 ±0.64 

0.069 0.75 ±0.10 

0.225 0.71 ±0.17 

0.841 0.86 ±0.50 

0.503 0.55 ±0.24 

0.162 0.72±0.15 

0.931 1.05±0.26 

0.935 1.07 ±0.36 

0.489 0.87±0.13 

0.819 0.96 ±0.25 

0.862 1.04 ±0.21 






,197 





.93 ± 


1 


.96 


0. 


.519 





.96 ± 


2 


,16 


0. 


.819 





.98 ± 





.11 





.474 





.94 ± 


2 


,75 



0.378 -0.25 ±0.29 

0.009 1.00 ±0.37 

0.682 -0.21 ±0.55 

0.629 -0.37±1.51 

0.009 0.90 ±0.36 

0.231 0.44 ±0.38 

0.981 0.01 ±0.42 

0.587 0.61 ±1.02 

0.746 -0.21 ±0.62 

0.036 0.49 ±0.24 

0.339 -0.40 ±0.41 

0.929 -0.07 ±0.68 

0.328 0.47 ±0.42 

0.603 0.18 ±0.34 

0.896 -0.05 ±0.42 

0.844 0.11 ±0.50 

0.985 -0.01 ±0.33 

0.055 -0.71 ±0.36 

0.696 0.13 ±0.33 

0.129 0.13±3.27 

0.454 0.07 ±3.58 

0.799 -0.04 ±0.16 

0.258 0.14±4.66 



0.436 0.85 ±0.05 

0.343 1.33 ±0.04 

0.201 0.65 ±0.23 

0.466 1.49 ±1.55 

0.720 1.12±0.08 

0.406 0.83 ±0.03 

0.543 1.20 ±0.55 

0.699 1.44 ±0.56 

0.199 0.48±0.24 

0.725 1.06±0.01 

0.447 0.81 ±0.06 

0.883 1.06 ±0.27 

0.407 0.48 ±0.29 

0.321 0.79 ±0.02 

0.046 2.04 ±0.09 

0.096 0.56±0.15 

0.507 1.18±0.02 

0.785 1.02 ±0.05 

0.532 0.87 ±0.04 

0.553 1.04±0.00 

0.546 1.04 ±0.00 

0.266 0.89±0.01 

0.088 1.16 ±0.01 



0.193 -0.18±0.15 
0.324 -0.18±0.19 
0.313 -0.16±0.16 
0.465 0.35 ±0.35 
0.800 0.04 ±0.18 
0.813 -0.04±0.16 

0.033 -0.27±0.13 
0.594 -0.15 ±0.25 

0.029 -0.40 ±0.18 

0.048 -0.19 ±0.09 
0.580 -0.08 ±0.15 
0.535 -0.11 ±0.16 
0.103 -0.24±0.12 

0.083 -0.23 ±0.13 
0.677 -0.06 ±0.14 

0.052 -0.45 ±0.19 
0.432 -0.10±0.12 
0.141 -0.14±0.12 
0.220 -0.15±0.12 

0.000 -0.14±0.03 
0.000 -0.13 ±0.04 
0.010 -0.15 ±0.06 
0.015 -0.11±0.05 



Note. - The fraction of absorption-dominated galaxies in each system is / a bs- A low value of P(t) 
indicates that two populations have significantly different means; a low value of P(F) indicates that they 
have significantly different variances. Sample-A consists of all systems except NRGsll7; Sample-D contains 
the systems with declining line-of-sight velocity dispersion profile, marked with a superscript "D" ; Sample-I 
consists of all systems except those in Sample-D and NRGsll7. 
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Table 5. Surface Density Profiles 



Spectral 7V gal NFW a 






NFW b 


Hernquist b 


Isothcrmal b 


Class X z / U 


c 




r c 






X*/v 




Sample-D c ; a p = 327-466 kms" 1 . 


Absorption 135 1.4 
Emission 75 2.1 
Combined 210 2.6 


4.9-12 

1.4-5.6 

4.3-9.1 


0.8 
0.60 
1.7 


0.062-0.16 
0.14-0.58 
0.074-0.16 


1.1 
0.79 
2.1 


0.26-0.46 

0.50-1.3 

0.30-0.49 


3.3 
0.79 
3.7 


8 Groups with a p > 350 km s 1 . 


Absorption 232 2.4 
Emission 131 1.5 
Combined 363 2.5 


7.6-16 

1.6-4.4 

4.3-7.7 


1.0 
1.0 
1.3 


0.067-0.14 
0.29-0.69 
0.12-0.22 


2.0 
1.1 
2.4 


0.30-0.46 

0.81-1.8 

0.45-0.65 


3.1 
0.9 
2.2 


9 Groups d with a p < 350 km s 1 . 


Absorption 102 5.7 
Emission 92 3.8 
Combined 192 4.2 


t 
t 
t 


4.5 
2.5 
3.0 


t 

0.17-0.59 
0.12-0.30 


3.7 
2.2 
2.5 


t 

0.52-1.53 
0.39-0.59 


8.8 
3.5 
5.1 



Note. — We list the 95.4% confidence interval for each fit parameter. N ga i is the total number of galaxies 
used in each fit. The core radius r c is always in units of h~ x Mpc. 

a Fit with projected radius R in units of r^o™ 1 . 

b Fit with projected radius R in units of Mpc. 

c The 5 Systems with Declining a p (R). 

d Includes only systems with > 10 members. 

'''This fit is ruled out with high confidence (x 2 jv > 3.5). 



